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Abstract 

We determine the maximal gap between the optimal values of an inte- 
ger program and its linear programming relaxation, where the matrix 
and cost function are fixed but the right hand side is unspecified. Our 
formula involves irreducible decomposition of monomial ideals. The 
gap can be computed in polynomial time when the dimension is fixed. 



1 Introduction 

We consider the general integer programming problem in standard form, 

Minimize c • u subject to A ■ u = b, u > 0, u integral. (1) 

Here A is a fixed d x n integer matrix, b e Z d and c G Q n . Its linear 
programming relaxation is obtained by dropping the integrality constraints, 

Minimize c • u subject to A ■ u = b, u > 0. (2) 

Suppose that the integer program (JTJ) is feasible and bounded. Then the 
linear program (J2J) is feasible and bounded as well, and the optimal value 
of (HJ) is greater than or equal to the optimal value of ©. We define the 
nonnegative rational number gap(v4, c) to be the maximum difference of 
the two optimal values as b ranges over vectors such that (0) is feasible 
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and bounded. It follows from known results [Tot Theorem 17.2] that this 
maximum is bounded. 

Our main aim in this paper is to provide an exact formula for gap(A, c). 
We express our results using the language of Grobner bases, as in |U Chapter 
8]) 0, [HI §4.4]. A nonnegative integer vector u = (ui,...,u n ) G N n is 
called non-optimal if it is not an optimal solution of ([TJ with b — A ■ u. 
We represent each non-optimal vector u by a monomial 
and we consider the ideal M(A, c) generated by these monomials in the 
polynomial ring . . . ,x n ]. The minimal generators of M(A,c) can be 

read off from a Grobner basis for If c is generic then M(A, c) is an initial 
ideal of the toric ideal of A; see [HI| • If c is not generic then we can compute 
M(A, c) using Algorithm 4.4.2]. A monomial ideal / in M[xi, . . . , x n ] is 
called irreducible if it is generated by powers of the variables: 

T / T u n+ 1 

1 \ -^ii ) x i 2 > • • • ' ^ir r 

Every monomial ideal M in M[xi, . . . , x n ] can be written uniquely as an 
irredundant intersection of finitely many irreducible monomial ideals, which 
are called the irreducible components of M. Suppose that / is an irreducible 
component of M(A,c). We define the gap value of / with respect to A and 
c to be the optimal objective function value of the auxiliary linear program 

Maximize w^c^ + u i2 c i2 H h u ir c ir — c-v 

subject to A ■ v = u h ■ A h + u i2 ■ A i2 H h u ir ■ A ir , (3) 

and v h ,...,v ir > 0. 

Here A\, A2, ■ ■ ■ , A n are the column vectors of the matrix A. 

Theorem 1.1. The integer programming gap, gap(A, c), equals the maximum 
gap value of any irreducible component I of the monomial ideal M(A,c). 

We remark that gap(A, c) is zero if and only if the monomial ideal 
M(A, c) is generated by squarefree monomials 

This paper is organized as follows. In Section 2 we rephrase our problem 
in the more general setting of lattice programs, and we prove Theorem II .11 in 
this context. In Section 3 we apply the work of Barvinok and Woods jSj on 
short rational generating functions to derive the following complexity result. 

Theorem 1.2. For fixed d and n, the integer programming gap, gap(v4,c) ; 
can be computed in polynomial time in the binary encoding of A and c. 
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Section 4 concerns applications to the statistical theory of multidimen- 
sional contingency tables. Here we are interested in the integer programming 
gap of certain higher-dimensional transportation problems. These play an 
important role in data security. For the statistical background see 0,0- 

In Section 5 we vary the cost vector c, and we prove that the function 

gap^ : R n ^R, c^gap(A,c) (4) 

is a piecewise-linear function on M n . We show that gap A is linear on the 
cones of a fan which refines the familiar Grobner fan (cf. jTTj and [T7]). 

We close the introduction with a small example. Let be the problem of 
making change using pennies, nickels, dimes and quarters, where the number 
of coins is fixed, and nickels and quarters are used most sparingly. In symbols, 



d = 2 , n = 4, A 



1111 
1 5 10 25 



c= [ 1 1 ] 



Using notation as in §1.5], this problem is solved by the Grobner basis 

g = { n 3 q - d 4 , nj - p 5 g , n 3 d* - p 5 q 2 , p 5 q 3 -d 8 }. 

The four underlined leading monomials generate the ideal 

M(A,c) = (n 3 ,p 5 ) n (n 3 ,q 3 ) H (n 6 ,d\q). (5) 

The gap values of the three irreducible components are 76/15, 4 and 5. Hence 

gap(A,c) = 76/15 = 5.0666666... 

This gap is attained when expressing one dollar and 14 cents with ten coins. 
The optimal solutions are (4, 2, 0, 4) and (0, 0, 136/15, 14/15) respectively. 

2 The gap theorem for lattice programs 

Lattice programs are defined as follows. Let £ be a fixed lattice of rank m 
in Z n . Then = £ ®i K is an m-dimensional vector space in W 1 . We also 
fix a cost vector c G Q n . Now, for every z G N n , we get a lattice program 

Minimize c • v subject to v = z mod £, v G N n , (6) 
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and its linear relaxation 



Minimize c • v subject to v = z mod £r, v > 0. (7) 

We define the lattice programming gap gap(£, c) to be maximum of the 
differences between the optimal values of (jOJ) and (|7jl as 2: runs over N n . 

Remark. The integer programs are lattice programs with £ = ker(j4)nZ n 
and 6 = A-z. The linear programming relaxations (J2J) correspond to the linear 
relaxations (0) of these lattice programs. Note that we have m = n — d. 

Let M(C, c) be the ideal generated by all monomials x u where u G N n is 
non-optimal for ©. Any irreducible component of M(C,c) has the form 

I(u,r) = {x^ ,x i2 2 ,...,x ir r ), 

where r = {ij, . . . , i r } and «6N" with Wj = for j ^ r. The gap value of 
I(u,t) with respect to £ and c is the optimal objective function value of 

Maximize u ix c ix + u i2 c i2 H h Uj r c ir - c-v . 

subject to v = u mod £k and u^, . . . , > 

We restate Theorem II. II for the more general setting of lattice programs. 

Theorem 2.1. The lattice programming gap, gap(£, c), equals the maximum 
gap value of any irreducible component I(u, r) of the monomial ideal M(C, c) . 

Example 2.2. (m = n) Let £ be a finite index sublattice of II 1 and c 
a nonnegative vector. Since £r = R n , the objective function value of ((7|) 
is always zero, so gap(£,c) is the largest objective function value of (@). 
The finite abelian group Z n /£ is in bijection with the set of monomials 
not in M(£,c). The irreducible components I(u, {1, . . . , n}) are indexed 
by monomials which are maximal with respect to divisibility 

among those not in M(C, c). The lattice programming gap is the maximum 
of the corresponding values c\Ui + C2U2 + • — h c n u n . □ 

Example 12.21 is called the group problem in the integer programming lit- 
erature P^l §24.2]. More generally, every lattice program © has a natural 
family of group relaxations which are indexed by subsets r of {1, 2, ... , n}\ 

Minimize c- v subject to v = z mod C, v G Z n , Vi > Wi G r. (9) 
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If an optimal solution v* of Q is a nonnegative vector then v* is also 
feasible and optimal for (JBJ). In this case we say that lattice program (0) 
is solved by r. The minimal collection of required subsets r is studied in 
[10J. The following result is well-known in the algebraic theory of integer 
programming; see §3]. The pairs (u,t) in Proposition 12.31 are called 
standard pairs. A combinatorial introduction can be found in [19]. 

Proposition 2.3. There is a unique minimal finite set S of irreducible ideals 
I(u, t) whose intersection is M(C, c) such that every optimal solution v* to 
(QJ) has the form v * = u + v ' , with v[ = for i G t, for some I(u, r) G S. 

When a lattice program © has an optimal solution v* = u + v' as in 
Proposition 12.31 we shall say that it is solved by the standard ideal I(u,t). 
In this case the group relaxation (jHJ) has the same optimal solution v*. 

Remark. The irreducible components of M(C, c) are a subset of S. In fact, 
the irreducible components are the minimal elements of S when the standard 
ideals are ordered with respect to inclusion. In the special case of Example 
l2~?l we have S = {l(u, {l,...,n}) : x u £ M(£,c)}. 

Lemma 2.4. Fix I(u,t) G S. The gap value of I(u,t) equals the maximum 
difference between the optimal values of (QJ) and Ffy as z ranges over all 
vectors in N n such that the program (0|) is solved by the standard ideal I(u, r). 

Proof. Suppose a lattice program © is solved by the standard ideal I(u, r) 
and has the optimal solution x* = u + x' where x\ = for all i G r. Let y* be 
an optimal solution for the linear relaxation (j7J). Then the difference between 
the optimal values is c ■ x* — c ■ y*. Since y* — x' is a feasible solution for (jSJ), 
the optimal value of this program is at least c-u — c- (y* —x') = c-x* —c-y*. 

Hence we only need to find a vector z G N™ whose associated lattice 
program © is solved by I{u, r), and such that the difference of the optimal 
values of (JUJ) and (JJJ) is greater than or equal to the optimal value of (JHJ. 
Let v* be an optimal solution to (jSJ) and define a vector v' G N n by v[ = 
max 

{°> - KM}- Then the lattice program JBJ with z = u + v ' is solved by 
I(u, r). In fact, z is the optimal solution, and v* +v' is a feasible solution for 
the linear relaxation (JJJ). The difference between the optimal solution values 
of the two programs is at least c • (u + v') — c • (v * + v') — c ■ u — c • v *. □ 

Proof of Theorem \2.1\ In light of LemmaEUand Proposition l2.3^ we just need 
to show that if I(u, r) and I(u', r) are two standard ideals with u < u', then 
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the optimal value of (JSJ) is at most that of (JHJ with u replaced by u' . In order to 
do this we will reformulate these programs. Let B be an n x m matrix whose 
columns form a lattice basis of C, and let {bi, . . . ,b n } C Z m be the rows of B. 
The feasible solutions to (JHJ) are in bijection with {t G M m : fcj-t < U{, i 6 r} 
via t i — > v = u—Bt. If we let w := C A 5 then the following linear program 

is equivalent to (jHJ) and has the same objective function value: 

Maximize w ■ t subject to bi ■ t < Ui Vi G r (10) 

The cost vector w is independent of u. If we replace u by u' in (jl(Jj) then the 
feasible region increases but the objective function is unchanged. Hence the 
optimal value of ()10j) can only increase when replacing « by «'. □ 

Example 2.5. Theorem 12. II suggests that we compute gap(£,c) by solving 
finitely many linear programs (jSJ), one for each irreducible component I(u, r) 
of M(C, c). One difficulty is that the number of irreducible components can 
be exponential in the problem size. We illustrate this phenomenon with an 
example taken from ^HJ §4] . For r > 4 let C r be the lattice generated by 

(r, r, r), (r — 1, r + 1, r — 1) and (0, 0, r — 2) in Z 3 . 

The index of C r in Z 3 is 2r(r — 2), so we are in the situation of Example 
12.21 Let c be a cost vector that chooses the degree lexicographically smallest 
solution of (jHJ). From the Grobner basis in fSJ Lemma 4.5] we find that 

M(A,c) = p| (x a , y\ z) H P (x, y\ z b ). 

l<a<r-2 l<a<r-3 
b=2r+l-a b=r-l-a 

The number of irreducible components is 2r — 5. This is exponential in 
0(log(r)), the bit complexity of the data. In the next section we demonstrate 
how the monomial ideals M(C, c) can be encoded more efficiently. □ 

3 Gap function as a short rational function 

In this section we prove Theorem II .21 The result extends easily to the lattice 
programs of Section 2, but for the sake of notational convenience, we present 
the relevant generating functions in the original setting of integer programs 
(JTJ). Throughout the section we assume that al and n are fixed integers. 
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Let NA denote the semigroup in Z d spanned by the columns A±, . . . , A n 
of the matrix A. The elements of NA are the feasible right-hand-side vectors 
of the integer programs we consider. Let flip : NA i— > Q be the function 
whose value at b is the optimal value of the integer program (0), and let 
fx lp : NA h> Q be the corresponding optimal value function for the linear 
program (J2J). Since c is assumed to have rational coordinates, we can com- 
pute (in polynomial time) a global denominator A 6 N by multiplying the 
least common multiple of the maximal minors of A with the least common 
denominator of c\, C2, . . . , c n . This choice of A ensures that every value of 
the functions \ijp or hlp is an integer multiple of A. 

We are interested in the generating function 

G(t x , ...,t d ;s) := J^tV^W-^W G Q[s 1/A ][[NA]]. 

beNA 

The ambient ring is the completion of the semigroup algebra of NA with 
coefficients in the univariate polynomial ring Qfs 1 ^]. We shall see that G is 
a rational series which can be represented by an element of 

Qfa, . . . ,t d )[s 1/A ] C Q(t u ...,t d ,s 1/A ). 

Thus G is a polynomial in s 1 ^ whose coefficients are rational functions in 
ti, . . . ,t d . The degree of that univariate polynomial is the integer program- 
ming gap g&p(A,c). We shall prove the following complexity result. 

Theorem 3.1. The rational function G(ti, . . . , t d ; s) can be computed in poly- 
nomial time in the binary encoding of the matrix A and the cost vector c. 

What we are claiming is that G(t; s) is a short rational generating function 
in the sense of Barvinok and Woods [3], and our proof is a direct applica- 
tion of their work. A different approach to the integer programming gap 
using generating functions and integer programming duality is presented by 
Lasserre in [13] . We illustrate the main point of Theorem 13 .11 in an example. 

Example 3.2. Let d — 1, n = 2, a, b £ N and consider the integer program 

minimize u\ subject to U\ + au 2 = b, u%,U2 > 0. 
The optimal values are Hlp(P) = and [iipib) = b mod a. Hence 

a— 1 00 a— 1 

G(t, s ) = j2t b s bmoda = J2J2 ti+jasi = X^/^ 1 ■ si - 

6eN 1=0 j=0 i=0 
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This is a polynomial in s with a terms. Its expansion requires exponential 
space in the bit size 0(log(a)) of the given data A = [1 a] and c = [10]. 
On the other hand, clearly this gap polynomial is a short rational function 

1 r .a+a, 

G(t,s) = 



(l-t a )(l - st)' 

Our strategy is to compute G(t, s) in polynomial time and then to extract 
gap(A,c) = degree s (G(t, s)) = a — 1. □ 

One ingredient in the proof of Theorem Kill is the Hadamard product 
* of two generating functions: if <7i(x) = J^Apc" and 52 (x) = '^2 v r )v'x v 
then <7i(x) * ^(x) is the generating function J2Pu7 u x u - The proof of [3J 
Lemma 3.4] tells us that, if <7i(x) and g<i (x) are short rational functions, then 
their Hadamard product pi(x) * <?a(x) can be computed in polynomial time. 
Another ingredient is the following lemma which is of independent interest. 

Lemma 3.3. The generating function for all the optimal points, 

H{x) = x " : u ^ optimal for (OJj with b = Au }. 

can be computed in polynomial time, in the binary encoding of A and c. 

Proof. The proof is an adaptation of §7.3]. Without loss of generality we 
assume that c G HI 1 . Let S = {u : x u G M(A, c)} be the set of non-optimal 
points, and let f(S;xi, . . . ,x n ) = J^ues^ ^ e ^ e generating function of S. 
This generating function is a rational function; in particular, we have 

n 

x) IJ(1 - *<) = g(S; x) 

i=l 

where ^(5*; x) is a polynomial. In view of the identity 

1 - g(S;x) 



H x 



nr=i(i-^) 



it suffices to show that g(S; x) can be computed in polynomial time. 

We let L := (n — d)D(A) where D(A) is the maximum of the absolute 
values of the maximal minors of A. Since we fix d and n, the bit size of L is a 
polynomial in the bit size of the data. Theorem 4.7 of PB] implies that if x m 
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is a term in g(S; x) then m 4 < L for % = 1, . . . , n. We let A := {(ui, . . . , u n ) G 
W 1 : u,i < L , V2} and /(A;x) its generating function. Furthermore let 

C := {(u,t;) G N 2n : Am = Aw, cu > cv + 1, and ttj < L, t> * < 2L, Vi} 

The projection of C onto the first n coordinates will be denoted by S', and 
we claim that 5" = S fl A. The inclusion S' C 5 fl A is clear. For the other 
inclusion let u G SnA. The theory of Grobner bases of toric ideals ^H] implies 
that there exists (u', v') G N 2n with u' < u, Au' = Av' and cu' > cv' + 1, and 
where u'-v^ = 0, u\ < L and v[ < L for i — 1, . . . , n. Now we let v = u — u' + v' 
and observe that (u, v) G C. Since S" is the projection of all lattice points in 
a polytope, Theorem 1.7 in [3] implies that /(S";x) = X!ues' x "> an< ^ hence 
g(S';x) = /(S";x) nr=i(l — x *)> can be computed in polynomial time. The 
claim we proved above says that 

/(5;x) = /(S';x) + ^ x", 

«gS\A 

and this implies g(S;x) = g(S';x) + fo(x), where h(x) is a series none of 
whose terms has its exponent vector in A. Now we conclude that 

g(S; x) = g(S;x) * /(A; x) = g(S'; x) * /(A; x). 

The Hadamard product on the right can be computed in polynomial time. □ 

Proof of Theorem \3.1\ We replace each variable X{ by t Ai s Ci in the rational 
function H(x). This monomial substitution can be done in polynomial time 
[3J Theorem 2.6]. The result is the short rational generating function 

H IP {t 1 ,...,t d ;s) = £tW 6 >. 

befiA 

The last series that is left to compute is 

H LP {t u ...,t d ;s) = J2 tbs ~" Lpib) > 

beNA 

since G(t; s) = i?/p(t; s) * t Hlp(^', s). Let cr = . . . ,i d } c {1, . . . , n} be 
an optimal basis of (0) for some b in NA The number of optimal bases a is 
constant since n and d are fixed. The union of S a := NA fl M>o{A : i G o~} 
as cr varies over all optimal bases is equal to NA. The generating function 
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f(S a ; t) = J2bes a ^ b can ^ e computed in polynomial time since S a is the 
set of lattice points in a rational polyhedral cone 0. We let c G R d be the 
unique vector such that {A 1 ■ c)j = q for i G a. Now the generating func- 
tion g ff (t;s) = Xlb6s CT ^ LP< ' 6 ' ) * s obtained from f(S a ;t) by the monomial 
substitution £j = tjS~ C4 for i = 1, . . . , d. Finally, we use [cO Corollary 3.7] to 
compute H LP (t] s) by patching the series g a (t; s) for the various bases a. □ 
We now show that gap(v4, c), which is the degree of G(t; s) as a polynomial 
in s, can be extracted in polynomial time. This uses the following lemma. 

Lemma 3.4. Let /(t; s) G Q(£i, . . . , id)[ s ] be a short non-zero rational func- 
tion and K a known upper bound on deg s (/(t; s)). Then this degree can be 
computed with log(-fT) Hadamard products o//(t; s) with polynomials in s. 

Proof. Without loss of generality we assume K = 2 k for some k G N. Let 



i=p 

and use the following binary search algorithm: 
f indDegree(/(t; s),p,r) 

0. If p = r return p. 

1. If /(t; s) * (7[ Pjr ](s) is not identically zero then 
f indDegree(/(t; s), |_(r + p)/2j , r). 

2. If /(t; s) * gip tr ](s) is identically zero then 
findDegree()(t;s), \p/2\,p-l). 

The call f indDegree(/(t; s), 0, i^) takes at most log(i^) steps to find the 
degree. Zero testing is done by substituting a positive vector (zi, . . . , Zd, w) 
where < z% <C 1 and < w C 1 which is not a pole of /(t; s) * g^ r ] (s) . 
Note that this Hadamard product is a polynomial in s whose coefficients 
have expansions of the form ^ 6gr t 6 for some set T. This means that the 
substitution of the above positive vector gives a positive number unless the 
Hadamard product is identically zero. □ 

Proof of Theorem \l .'A By Theorem 13 .11 we can compute the rational function 
G(t; s) in polynomial time. The degree of G(t; s) viewed as a polynomial in s 
is gap(v4,c). Theorem 17.2 of [T3j implies that gap(A,c) < nD(A) XT=i l c *l- 
Let K be this upper bound. We observe that log(f^) is a polynomial in the 
bit size of the data. We can hence use Lemma [3.41 to find gap(v4, c). □ 
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4 An Application to Algebraic Statistics 



We present an application to the statistical theory of disclosure limitation. 
See jS] and P and the references therein. Suppose we are given data in 
the form of an n- dimensional table of nonnegative integers. The aim is to 
release some marginals of the table but not the table's entries themselves. 
If the range of possible values that a particular entry can attain in any 
table satisfying the released marginals is too narrow, or even worse, consists 
of the unique value of that entry in the actual table, then this entry may 
be exposed. This shows the importance of determining tight integer upper 
and lower bounds for each entry in a given table. A choice of marginals 
corresponds to a simplicial complex on {1,2, ... ,n} and can be represented 
by a zero-one matrix A, as described in jHJ §1]. In statistical language, the 
matrix A specifies a hierarchical model for a contingency table with n factors. 

The table entry security problem can be formally stated as follows: sup- 
pose u is a table with nonnegative integer entries, where the marginals are 
computed according to a fixed hierarchical model A and let u^...^ be a par- 
ticular cell for such tables. Compute optimal lower and upper bounds L and 
U such that L < u^...^ < U for all tables with the same marginals as u. 

This problem is an integer program ([!}: minimize (or maximize) m^^...^ 
over all tables with nonnegative integer entries subject to fixing the marginals. 
In view of the difficulty of solving integer programs exactly, various re- 
searchers resorted to solving the linear programming relaxation (J2J) instead: 
minimize (or maximize) Ui 1 i 2 ...i n over all tables with nonnegative real entries 
subject to fixing the marginals. This relaxation is tractable, but it usually 
fails to deliver the exact integers L and U. One faces the problem of esti- 
mating the integer programming gap for the table security problem. 

We give the precise definition of the relevant matrices A. Consider d\ x 
■ ■ • x d n -tables with entries u^...^ where 1 < % j < dj. We fix a hierarchical 
model by specifying a collection of subsets Fi, . . . , of {1, . . . , n}. The 
marginals of our table are computed with respect to these subsets. If Fi = 
{ii-, ■ ■ ■ ijs} then the Fj-marginal is a x • ■ • x dj s table b with entries 

b kl -k s = Yl u h-in- ( n ) 

Example 4.1. The classical transportation problem corresponds to di x ea- 
tables where the marginals are computed with respect to F\ = {1} and F 2 = 
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{2}. The three-dimensional transportation problem concerns dx x d 2 x ea- 
tables with Fi = {1, 2}, F 2 = {1, 3}, and F 3 = {2, 3}. The marginals are 

6 ifc =x;^. 

fc i i 

For a discussion from the Grobner basis perspective see ^01 §14. C]. □ 

We define A to be the zero-one matrix with d\d2---d n columns that 
corresponds to the linear map that computes the marginals of tables. We 
let u be the vector of variables representing the cell entries. Then A ■ u is 
the vector of the k lower-dimensional tables computed as in There is 

a transitive symmetry group acting on the columns of A, so it suffices to 
examine the particular cell entry Uxx-i which corresponds to the first column 
of A. The table entry security problem is the pair of integer programs 

Minimize (Maximize) Uxx-i subject to A-u = b, u > 0, u integral. (12) 

Our Theorem 11.11 gives an exact formula for the integer programming gap 
of these problems, which we abbreviate by gap_(A) and gap + (A) respectively. 
Thus gap + (A) is the worst error one gets when using linear programming in 
computing the bound U for any d\ x • • • x d n -table with any fixed margins b. 

We illustrate our results for 2 x 2 x 2 x 2-tables (wy-jw). The K^-model is 
specified by taking all six two-dimensional margins Fx — {1, 2}, F 2 = {1, 3}, 
F 3 = {1,4}, F 4 = {2,3}, F 5 = {2,4}, F 6 = {3,4}. The zero-one matrix A 
for the -fQ-model has 24 rows and 16 columns: 



A = 
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Here the cell entries are ordered lexicographically, Wim, W1112, • • • ,^2222- 
We have the following computational result. 

Proposition 4.2. Every hierarchical model for 2 x 2 x 2 x 2-tables satisfies 

gap + (A) < 5/3 = 1.666666..., 
and this bound is attained for the K^-model, that is, for the matrix A above. 

Proof. We show that gap + (A) = 5/3 for the .fQ-model. Similar computations 
show that gap + (y4) < 5/3 for all other simplicial complexes on {1,2,3,4}. 

The monomial ideal M(A, c) where A is the above matrix and c = 
[— 1, 0, . . . , 0] has 61 minimal generators. Two of them are 

a '1112 a '122l3 ; 12223 ; 212l3 ; 21223 ; 221l2 ; 2212 and ^1122^1212^1221^2111^2222- 

This ideal has 139 irreducible components. One of these components is 

I(u,t) = (Xin2) ^1121) ^1122) ^1211) ^1212) ^1221; ^1222; 

2 2 \ 

X 211lJ X 2112) ^2121? ^2122) ^2211 j ^2212 j ^2221; ^2222/- 

Here the non-zero components of u are Wim = M1121 = M1211 = M2111 = 
^2222 = 1- The gap value of I(u, r) is 5/3. We also give the optimal solution 
v for the linear program (j2J) with b = A ■ u: 

vni2 = v n2 i = v 12 n = f 2111 = ^2222 = 0, 

"^1122 = "^1212 = "^1221 = "^1222 = ^2112 = 1/3, 
^2121 = ^2122 = ^2211 = ^2212 = ^2221 = 1/3, V im = 5/3. 

□ 

The following corollary was pointed out to us by Rekha Thomas. 

Corollary 4.3. Let M(A, c) be the monomial ideal corresponding to the min- 
imization problem i fTHj) . Then gap„(A, c) + 1 is at least the maximum degree 
of Xn... i in any minimal generator of M(A,c). 
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Proof. This statement is equivalent to 

gap_(A) > max{%...i : I(u,t) irreducible component of M(A,c)}. 

Let I(u,t) be an irreducible component such that U\\...\ > 1. Since c = 
[1, 0, . . . , 0], the objective function of the program (|TU|) corresponding to this 
component would be bu...yt. Moreover, by Theorem 2.5 in JU] 5 the inequality 
&11-1 ■ t < Wii-i is a facet of the feasible region of (flDjl. Hence the optimal 
solution value is Mu.„i, Now Theorem 11.11 gives the result. □ 

Remark. In the statistics literature there are various approaches to estimate 
L and U of the table security problem. For general hierarchical models an 
iterative algorithm for such an estimation is given in [I]. A detailed analysis 
for decomposable models is given by Dobra and Fienberg 

5 The gap fan 

In this section we allow the cost function c to vary in the programs or 
(jSJ). For each fixed matrix A (resp. fixed lattice C) we thus get a function 

gap A : W n — ► R, c i — ► gap(A, c). (13) 

Our goal is to show that this function is piecewise linear, and the natural 
pieces on which the function is linear form a fan which we call the gap fan. 

Consider the Grobner fan of the matrix A. Following jTHj and ^7], this 
is the coarsest polyhedral fan in M n on which the map c i— > M(A, c) is 
constant. Efficient software packages for computing the Grobner fan are de- 
scribed in JTI] and ^2j- The gap fan will be a refinement of the Grobner fan. 
Hence it suffices to describe the gap fan on each Grobner cone K separately. 

We fix a maximal cone K, of the Grobner fan. By the results of J7j, the 
polyhedral cone /C consists of cost vectors such that the optimal solutions of 
all integer programs (0) are constant as the right-hand-side vector b varies. 
There exists a monomial ideal M in R[xi, . . . , x n ] such that 

K, = {ceM. n : M(A,c) = M}. 

Let {l(u, r)| be the finite set of all irreducible components of the monomial 
ideal M. For each such component let v UjT be the optimal solution to the 
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linear program (j3J) . It is not hard to see that these optimal solutions are also 
constant as c varies in JC. Now we define the following polyhedral cone: 

Q ■= [( c ,t) G R n+1 : c G K and c- (w - t> MjT ) < t V J(u,r)} (14) 

Theorem 5.1. JTie projection of the lower envelope of Q onto K, is a poly- 
hedral subdivision of K,. The gap function is linear on each face of this 
subdivision. 

Proof. The facets of the lower envelope of Q are of the form {(c, t) G M n+1 : 
c • (u — v u>T ) = t} where t is the optimal value of the program Q for all c 
such that (c, t) is on this facet. For such pairs (c, t) we have t = gap(A, c), by 
Theorem ll.il Hence the lower envelope of Q is the graph of the gap function 
(Unj) over /C. Now, clearly, the projection of the lower envelope onto /C is a 
polyhedral subdivision of /C, and by its construction, the gap function (|T3|) 
is linear on each cone in this subdivision of /C. □ 

We define the gap fan as the refinement of the Grobner fan of A, where 
each Grobner cone /C is subdivided as in Theorem l5.ll Our discussion implies: 

Corollary 5.2. The gap function Mty is a piecewise linear function on M n . 
It is linear on the cones of the gap fan, but generally not on the Grobner fan. 

Example 5.3. Let us return to the example in the Introduction but now with 
varying cost function. Our problem is to make change using pennies, nickels, 
dimes and quarters, where the number of coins is fixed. The Grobner fan of 
the corresponding 2 x 4-matrix A has seven maximal cones. These cones and 
the corresponding ideals can be computed using TiGERS [11 or CaTS [12J. 

The gap fan has eight cones. Exactly one cone of the Grobner fan is 
divided into two cones. It is the one corresponding to the ideal in ©. This 
Grobner cone is defined by the inequalities 

3n + q > Ad and 8d < 5p + 3q. 

The hyperplane defined by 

305p - 135n - 308d + 138g = 

splits this cone into two pieces. On the positive side of the hyperplane the 
gap(v4, c) is given by the irreducible component (p 5 , n 3 ), and on the negative 
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side it is given by (n 6 , d 4 , q). We list the irreducible components of all seven 
initial ideals and the winning irreducible component for each of them: 



M(A,c) 


Winning component (s) 


(p 5 ,d 4 ) 


(p 5 ,rf 4 ) 


(p 5 ,q) H (p 5 ,n 3 ) H (d 4 ,q) 


(p 5 , n 3 ) 


(p 5 , n 3 ) [~l (n 9 , q) 


(p 5 , n 3 ) 


(p 5 , n 3 ) n (n e , q) H (n 3 , q 2 ) 


(p 5 , n 3 ) 


(p b , n 3 ) n (n 6 , d A , q) H (n 3 , g 3 ) 


(p 5 , n 3 ) 
(n 6 ,d\q) 


{n\d\q) n (n 3 ,d 8 ) 


(n 6 ,d\q) 


(n 6 , rf 4 ) 


(n 6 , d A ) 
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